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Abstract: We consider the most general three-state spin chain with U(l) 3 symmetry 
and nearest neighbour interaction. Our model contains as a special case the spin chain 
describing the holomorphic three scalar sector of the three parameter complex deforma- 
tion of Af = 4 SYM, dual to type IIB string theory in the generalized Lunin-Maldacena 
backgrounds discovered by Frolov. We formulate the coordinate space Bethe ansatz, calcu- 
late the S'-matrix and determine for which choices of parameters the S'-matrix fulfills the 
Yang-Baxter equations. For these choices of parameters we furthermore write down the 
i?-matrix. We find in total four classes of integrable models. In particular, each already 
known model of the above type is nothing but one in a family of such models. 
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1. Introduction 

Integrability has played a prominent role in recent years exploration of the AdS / CFT cor- 
respondence [0], the starting signal being the discovery by Minahan and Zarembo that the 
planar one-loop dilatation operator of the so (6) sector of N = 4 SYM could be identified 
with the Hamiltonian of an integrable spin chain M. This discovery was soon extended 
beyond the so (6) sector || and beyond the one-loop order §. On the string theory side inte- 
grable structures first made their appearance with the observation that the Green Schwarz 
superstring on AdS^ x S 5 possessed an infinite set of non-local conserved charges || and 
were more precisely described in a series of important papers || |?], ^, ^, [H], [llfl . An inter- 
esting question is to what extent the existence of integrable structures in the gauge-string 
duality is linked to super symmetry. 1 A first possibility to study this question was provided 
with the discovery by Lunin and Maldacena of the background required to make type 
IIB string theory dual to the Leigh-Strassler /3-deformation of Af = 4 SYM [jl4| , a confor- 
mal quantum field theory carrying M = 1 supersymmetry. Later the construction of Lunin 
and Maldacena was generalized to a three parameter complex deformation of J\f = 4 SYM 
which is still conformal (at least to leading order in N) but for which supersymmetry is 
lost [[Tq ]. Furthermore, it was shown that for uni-modular deformation parameters one can 
construct a Lax pair for (the bosonic part of) the corresponding classical string theory, 
showing that this theory like its undeformed counterpart is integrable pl| . (A Lax pair for 

1 It is of course also interesting to study the existence of integrable structures for gauge theories for which 
no string theory dual is currently known. For a review and further references, see [jL2|. 
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the Landau-Lifshitz model describing fast rotating three-spin strings in the deformed back- 
ground [16] was derived in reference fl7|| ). Let us denote the three deformation parameters 



as rje 1 ^, with rj > 0, jj real, j € {0, 1,2}. (For the precise definition, see section ||.) The 
three parameter deformation of N = 4 SYM reduces to the Leigh-Strassler one when the 
three deformation parameters are identical. In [18| the general deformation process was 



studied from the gauge theory point of view and it was found that the case vq = r\ = T2 = 1 
always leads to an integrable dilatation operator (at least to the order that the latter is 
known). In addition, exact matching of gauge theory anomalous dimensions and string 
state energies was found in the two- as well as the three-scalar holomorphic sub-sectors 



(the analogues of the su(2) and su(3) sub-sectors of the undeformed theory) |fl6 |. Results 
on the existence of integrability on the gauge theory side for general deformation param- 
eters are mostly negative. For instance, in the two scalar holomorphic sector integrability 
holds at one loop order for ro = r\ = r2 = r ^ 1 and 7o = 7i = 72 = 7 but is very likely to 
break down already at two-loop order |fl9|| . Furthermore, it was found in [20|, see also plf , 



that assuming ro = r\ = r2 = r integrability can only be obtained at one loop order in 
the three scalar holomorphic sector if r = 1. The one- loop dilatation operator in the three 
scalar holomorphic sector of the three parameter complex deformation of the gauge theory 
is a three-state spin chain with nearest neighbour interaction and U(l) 3 symmetry. In the 
present paper we study the most general spin chain with these properties and determine 
for which choices of parameters it can be integrable. We find four classes of integrable 
models. In particular each already known integrable model of the above type is nothing 
but one in a family of such models. We also rule out the possibility that the complex three 
parameter deformation of the one-loop planar su(3) sector of Af = 4 could be integrable 
under more general circumstances than those already identified. 

The organization of our paper is as follows. We begin in section ^ by writing down 
the most general Hamiltonian for the type of model we wish to consider. Subsequently, in 
section ^ we determine its S-matrix and in section |] investigate for which choices of param- 
eters this S-matrix fulfills the Yang-Baxter relation. For these choices of parameters we 
then in section [| construct an i?-matrix for the model. Finally, we compare our integrable 
models to those already known in section |(| We finish with a brief conclusion. 

2. The model 

We consider the most general 3-state spin chain with nearest neighbour interaction and 
U(l) 3 symmetry. We denote the basis of states at a given site as |0), |1), and |2). A U(l) 3 
symmetric Hamiltonian conserves the number of states of each type and accordingly takes 
the form 

H = HqqEqoEqo + H\\EuEu + H22E22E22 + HI2E11E22 + #21 ^12-^21 
+Hl 2 V E2\Ei2 + H21E22E11 + hIiEiqE^i + hIqEuEqq + H^\EqqEh 
+H®qEqiEio + HQ2E20E02 + H2QE22EQ0 + Hq 2 2 EooE22 + H2QE02E20, (2.1) 

where E^E\ m is an abbreviation for ^2f = i E^E 1 ^ 1 , with L being the length of the spin 
chain, and where Eij\k) = \i)5jk, i-e. are the generators of gl(3). The spin chain is 
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assumed to be closed. Restricting ourselves to Hermitian Hamiltonians we shall assume 
that the diagonal elements are real and that the off-diagonal ones can be written as 



Hll = {HUT = r e^°, H™ = « 2 )* = ne^ , H§ = « )* = r 2 e^ , (2.2) 

where ro, r\ and r 2 are real and positive and where 70, 71 and 72 are real. Our aim 
is to determine for which values of the parameters this Hamiltonian is integrable. An 
interesting spin chain which is included in the class of models given by eqn. (|2.1| ) is the 
spin chain describing the the three scalar holomorphic sub-sector of the three parameter 
complex deformation of M = 4 SYM, dual to type IIB string theory on the generalized 
Lunin-Maldacena backgrounds found by Frolov. For this model one has 

trOO _ rrll _ tj-22 _ n rrl2 _ rjOl _ tj-20 _ -, 

rr21 _ r 2 tt02 _ 2 ttW _ 2 t<y o\ 

n 2\ — '0' : -"02 — '1' -"10 — '2' 1^-°/ 

H\l = (Hll)* = r e**> , H§ = « )* = ne^ , H™ = « 2 )* = r 2 e^ . 



In [20 1 it was found that is integrable for 70 = 7i = 72 = 7 = 2 t zl with n integer 



and ro = T\ = r 2 = r only if r = 1. Later the model was shown to be integrable for 



fQ = T\ = r 2 = 1 for arbitrary values of 70, 71 and 72 [18]. The case = r\ = r 2 = 1, 
7o = 7i = 72 = 7T is the usual XXX su(3) spin chain. An additional number of integrable 
special cases of the model (2.1) are already known. We shall return to these in section |6|. 



The problem of studying the most general spin chain with certain global conservation laws 
has also been pursued in condensed matter physics, see for instance p4{] . Very recently, 
another problem of a somewhat similar nature was addressed, namely the problem of 
determining the most general integrable long range spin chain with the spins transforming 
in the fundamental of fll(n) ]pq] , 

Integrability implies the existence of an i?-matrix, R(u), depending on a spectral pa- 
rameter u and fulfilling 



R(u)\ u=U0 = V, V—R{u)\ u=U0 = H, (2.4) 



d_ 

du 

where V is the permutation operator. Furthermore, the i?-matrix must obey the Yang 
Baxter equation 

Rj^ia r)/^;(,,;/^;ir] /^(/•l/^^.l/^i^ r). (2.5) 



We can write a Hamiltonian H of the form ( |2.1| ) with the angles 70, 71 and 72 present in 
terms of a Hamiltonian without angles, H, 

Hfj = exp Q(ey m 7 m - e fc/n 7 n )^ H§. (2.6) 

Then, if the Hamiltonian H admits an .R-matrix, Rfj, we can construct an i?-matrix, R^j, 
corresponding to H by the following recipe [26, 18 1 



R\\ = exp ( ^(e iim 7 m - e fc/n 7„) ) R l £. (2.7) 
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Conversely, given an i?-matrix corresponding to a Hamiltonian of the form ( |2,l[) with the 
angles set to zero, one can construct an i?-matrix for a Hamiltonian with the angles re- 



introduced by applying the transformation inverse to (2.7). Thus any model integrable 
without the presence of angles is also integrable with arbitrary values of the angles and 
any model integrable with non-zero values of the angles is also integrable if the angles are 
set to zero. In other words, when searching for integrable models, one can leave out the 
angles from the analysis. 

To test for integrability in a system that admits scattering one can, instead of directly 
constructing an /2-matrix, study the properties of an appropriately defined S-matrix. In 
the present case the S-matrix can be written down by choosing a reference state, say, 
|0. . . 0} and considering scattering of excitations of type 1 and 2. Integrability normally 
requires non-diffr active and factorized scattering. 2 This implies that the S'-matrix must 
fulfill the Yang-Baxter like relation 

S2,3{P2,P3)Sl,3{pi,P3)S 1:2 (pi,P2) = Sl j2 (pi,P2)Sl > 3(pi,P3)S2,3{P2,P3), (2-8) 

where each S'-matrix acts in a eight-dimensional space. For a graphical representation of 
this relation, see Appendix A. In the following we shall follow the simpler route of first 
studying the S-matrix. Hence we start by considering two particle scattering, write down 
the coordinate space Bethe ansatz, determine the S-matrix and investigate for which choices 
of parameters it fulfills the relation fl2.8| ). Afterwards, we write down the corresponding 
i?-matrix. The coordinate space Bethe ansatz was introduced by Bethe in 1931 123] and 



revived by Staudacher in connection with the study of J\f = 4 SYM in 23]. In reference [24] 
a model similar to ours was analyzed using the so-called matrix product ansatz, but the 
analysis was not taken to the point of actually determining and classifying all integrable 
cases. 



3. The Bethe Ansatz and the S-matrix 

It is obvious that the three states |0 . . . 0), |1 . . . 1) and |2 . . . 2) are all eigenstates of the 
Hamiltonian. In the following we shall choose |0 . . . 0) as our reference state. The integra- 
bility properties of the model should not depend on the choice of reference state. First, let 
us study states containing only one excitation. We define 

|1) = Yl i 00 °)' 

l<h<L 

and similarly for |2). Writing down the Schrodinger equation gives 

(H™{L - 2) + H$ + H°l) Mh) + n (Mh + 1) + Mh ~ 1)) = ElMh)- (3-2) 
This equation is immediately solved by the plane wave 

Mh) = e ip ^, Ei = H ™(L - 2) + H™ + fl™ + r x (e ipi + e" ipi ) . (3.3) 

2 In the case where the angles are present in ( |2,2| ) the scattering is actually diffractive, see appendix A. 
Due to the argument given above, the model can nevertheless be integrable. 
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Similarly, we find a one excitation eigenstate |2) given by 



JP2h 



E 2 = H° ° (L - 2) + H™ + H° 2 2 + r 2 (e^ + e"* 2 ) 



(3.4) 



We notice that in general the two types of excitations have different dispersion relations. 
Whereas a difference in the p- independent terms is harmless, a difference in the p dependent 
terms is normally viewed as a signal that the model is not integrable. We shall make this 
statement more precise below. 

Let us now turn to studying two-body interactions. We choose as a basis for the space 
of two-particle excitations the states |11), |12), |21) and 1 22) where \ij) is defined as 



Aj(h, h)\00 i 000 j 000). 



(3.5) 



l<h<h<L 



When ro / the model allows for scattering between particles of type 1 and 2 and the 
states 1 12) and 1 21) will mix. Accordingly, our S-matrix will take the form 



S 



(a 0\ 

c b 

b c 

\0 dj 



(3.6) 



We shall first consider the case ro / 0. Later, we will comment on the simpler case ro = 0. 

Let us to begin with study the states |11) and 1 22) which do not mix with anything. 
For the former one we get from the Schrodinger equation 



h > h + 1 : 



«°(L - 4) + 2H™ + 2< 1 ) ^{hM) 

+n frMZi + i,z 2 ) + V'ii(M2 + i) 

+^ll(il-M2)+^ll(*l>*2-l)}, 



(3.7) 



h = h + l: 

En^lihM) = (H$(L - 3) + H\\ + H§ + H*\) ^{hM) 
+n frMZi, Z 2 + l) + (ii - l, fe)} • 

Inserting the standard Bethe ansatz 

^n(ii,i 2 ) = e ipiil+ ^ 2 + afe,^)^ 1 ' 2 ^ 1 , 

we find 

E n = H™{L - 4) + 2^° + 2# ° 1 1 + n (e* 1 + e ip2 + e~ ipi + e^) 

and 



a(Pi,P2) = - 



aie i P i + rie »pi+«P2 + ri 
aie iP2 + rie ipi+ip2 + ri ' 



(3.8) 



(3.9) 

(3.10) 
(3.11) 
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where 

ai = - H$ - H\{ + Hq\. (3.12) 

The results from 22-scattering can be obtained from those of 11-scattering by making the 
replacements H\q — > H 2 ®, HqI — ► H® 2 , r\ — > r 2 . Thus we use the Bethe ansatz 

il>22(hM) = e^ h+ip2h + d(p 2 , Pl )e ip ^ +ip2h , (3.13) 

and get 

E 22 = H™{L - 4) + 2H™ + 2tf ° 2 2 + r 2 (e ipi + e^ 2 + e^ Pl + e"^ 2 ) , (3.14) 

and 
where 

°2 = #20 _ #00 _ H M + #02 • (3.16) 

Next, let us turn to studying the states which mix. For these the Schrodinger equation 
gives rise to the relations 

h > h + 1 : 

E^ 12 {hM) = {H$(L - 4) + H™ + #o°i + # 2 2 ° + #o°f ) ^12 (/i, fe) (3.17) 
+ri (Mh + 1, fe) + Mh - 1, fe)) 

+r 2 ($1 2 (*1, + 1) + ^12(^1,^2 - 1)) , 

E^ 21 (h,l 2 ) = (H™{L-A) + HlS + H*l + H™ + H™)i> 21 (hM) (3.18) 
+n ($21 (ii, Z 2 + 1) + $>i(/i, Z 2 - 1)) 

+r 2 ($21 (ii + 1, i 2 ) + $21 (*1 " 1, fe)) , 

= ii + 1 : 

EiP 12 (h,l 2 ) = {H™(L - 3) + ffJ! + fl™ + F| °) $ 12 (Z 1 , / 2 ) (3.19) 
+n>$2i(/i,/2) + ri$i 2 (Zi - l,/ 2 ) + r 2 ipi 2 (h,l 2 + 1), 

EiP 21 (h,l 2 ) = (F ° °(L - 3) + + + # °f) $2i(ii, Z 2 ) (3.20) 
+ro$i2(/i,/2) + riip 2 i(h,l 2 + 1) + r 2 ip 21 (h - l,l 2 ). 

Let us now choose the following Bethe ansatz 

^12 = ,4 12 e^ lZl+ip2Z2 + i4i 2 e ip, i ,2+i ^ ,1 ) (3.21) 
$21 = ^ 21 e ip 'i /l+ip 2 /2 + i4' 21 e ipi,2+ipa,1 ) (3.22) 

where due to the translational invariance of our model 

Pi+P2=p[+p' 2 - (3.23) 
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The idea is that two excitations with momenta p\ and p 2 can scatter whereby their momenta 
get changed to p\ and p 2 . Inserting this into ( [3. 17 ) and ( 3.18] ) we find for the energy 



E = H™(L - 4) + ffJo + H$ + + fljg + rife* 1 + e" ipi ) + r 2 (e^ + e"^) 
= Hqq(L - 4) + + # i + F| ° + tf °2 2 + n(e^ + e"^) + r 2 (e^ + e~^X3.24) 

Equations ( |3.24| ) and ( |3.23| ) determine and p 2 

e %p i=e tpi 2 1 . ^. , (3.25) 
e ip 2 = e *w lULlgl , (3.26) 

where the solutions are chosen such as to reproduce the standard case p[ = p\, p 2 = P2 
when r\ = r 2 ^ 0. If both r\ and r 2 vanish the usual Bethe ansatz p[ = p\ and p' 2 = P2 
is still applicable. If only one of the two vanishes one necessarily has p' x = pi and p' 2 = pi, 
i.e. scattering between excitations of type 1 and 2 is not possible. In this case one should 
choose another reference state. However, the results for this situation can be obtained by 
symmetry arguments from those of the case ro = which we will consider below. 

The S-matrix elements involved in the 12-scattering are defined, using the transmission 
diagonal representation, as 

A 2i \ = ( c(p2,Pi) b(p 2 ,Pi)\ ( A X2 \ 

A' l2 ) \b( P 2,Pi) c(p 2 , Pl )J \A 21 J> [ ■ ' 



and can be found from ( gl9|) and (fjg) which with the Bethe ansatz (|372l|) and (|3T2^ ) 
read 

= A 21 r e ip * + A' 21 r e ipi (3.28) 
-A 12 {n + r ie ipi + r 2 e~ ip2 } e ip2 - A' 12 {n + ne ip i + r 2 e~ ip ^ } e ip 'i , 

= A 12 r e ip2 + A' 12 r e ip 'i (3.29) 
-A 21 {r 2 + ne"^ + r 2 e 4p 'i} - ^ {r 2 + ne"* 1 + r 2 e* P2 } e ip \ 

where 

r i = #io ~ #oo - #12 + #02 > (3.30) 

72 = #°i - #oo - #21 + #20- (3.31) 

As mentioned above it is common lore that the model can not be integrable unless the two 
scattering excitations have the same dispersion relation. We will give the precise argument 
for our Hamiltonian in Appendix A. 3 ^From now on we assume that r\ = r 2 . If in addition 
rp 7^ scattering will also be possible with the choice of either of the states |1 . . . 1) or 



3 Integrable Hamiltonians of the type ( |2.l[ ) with generic values of the angles, however, constitute an 
exception to the rule, see appendix A. 
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1 2 . . . 2) as reference state and it follows by symmetry arguments that ro = r\ = r 2 = r. 
For ro = r\ = r 2 = r ^ the remaining S-matrix elements read 

c (p 1}P2 ) = c( Pl ,p 2 ) = 1 ( e *w - (1 + e^+^) , (3.32) 

b(pi,J>2) = -^((*i*3 " l)e ipi+ip2 + + 

+ (tie ipa + i 2 e ipi ) (1 + e *i-HB»)) , (3.33) 

Kpi,P2) = ~^((tih - l)e^+^ + (l + e *i-N») 2 

+ (i 2 e 4p2 + he ipi ) (1 + e^^)) , (3.34) 

where 

D(pi, P2 ) = ((ht 2 - l)e 2ip2 + (1 + e ipi+ip2 f + (ti + i 2 )e* P2 (l + e ^+^)) , (3.35) 

and 

*i = n/r, t 2 = r 2 /r. (3.36) 
It is easy to verify that the S-matrix given by the relations ( |3.11|) , (|3.15| ), ( |3.32|) , (|3.33| ) 



and ( 3.34[ ) is unitary. We notice that effectively the S-matrix (with the restriction ro 



T\ = r 2 = r 7^ 0) depends only on the four parameters 

si = a 1 /r, s 2 = a 2 /r, ti, t 2 , (3.37) 



(cf. eqns. ( gig ), ( fUBD , ( p2|) , ( g3g ), ($$M) and (p5|)). This is easy to explain. To begin 



with we had 15 parameters, we then removed the angles as integrability properties can be 
analysed without them and set ro = r\ = r 2 = r since the model otherwise is necessarily 
non-integrable. This leaves us with 10 parameters, the single off-diagonal one, r, and 9 
diagonal ones. First, one can of course make a global rescaling by 1/r without changing 
the scattering properties of the model. Secondly, one can construct the following number 
operators 

(3.38) 
(3.39) 
(3.40) 

where iVj and Mi counts the number of particles of type i. Only five of these operators are 
independent since we have the relation 

N + iVi + N 2 = M + Mi + M 2 . (3.41) 

Since the number operators commute with the Hamiltonian, adding such operators will not 
change the scattering properties of the system. This means that effectively the S-matrix 
depends only on four parameters. 



N 


= 1(5 


5 -Ebo> 


M 


= E 00 $ 


5 1 


Ni 


= li 




Mi 


= E n a 


5 1 


N 2 


= 1<S 


5 E 22 , 


M 2 


= E 22 $ 


5 1 
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When n = r 2 = we have 



c{pi,p 2 ) = c(pi,p 2 ) = b(pi,p 2 ) = b(p 1 ,p 2 ) = a(pi,p 2 ) = d(p 1 ,p 2 ) = -e tpi W2 . (3.42) 
Furthermore, for ro = 0, r\ = r 2 = r the S-matrix elements involved in 12-scattering 

read 

c(pi,p 2 ) = c(pi,p 2 ) = 0, (3.43) 
b(p uP2 ) = - 2 . —. -. + . (3.44) 

b(pi,p 2 ) = — ^ ^ • (3.45) 

4. Solution of the YBE's 

A necessary condition for integrability is that the unitary S-matrix fulfills the Yang-Baxter 
equation 

S 2 ,3(P2,P3)Sl,3{pi,P3)S lj2 (p 1 ,p 2 ) = S lj2 (p 1 ,p 2 )S l! 3(p 1 ,P3)S 2j 3(p 2 ,P3), (4.1) 

where each S-matrix acts in a 8-dimensional Hilbert space. 

For a S-matrix on the form ( |3.6| ) the Yang-Baxter equation gives rise to 7 relations 
between matrix elements, 

bb'b" - bb'b" = 0, (4.2) 

ab'a" - ba'b" - cb'c" = 0, (4.3) 

ca'b" - ac'b" + bb'c" = 0, (4.4) 

be a" - cb'b" - ba'c" = 0, (4.5) 

db'd" - Mb" - cb'c' = 0, (4.6) 

cd'b" - dc'b" + bb'c" = 0, (4.7) 

be' d" - cb'b" - bd'c" = 0, (4.8) 

where b = b(pi,p 2 ), b' = 6(^1,^3), b" = b{p 2 ,p^), etc. At first sight these equations seem 
rather involved but a systematic investigation with the aid of Mathematica allows one 
to find the most general solution by purely analytical means. The simplest equation is 
equation ( |4.2j ) so this is the one to address first. Next, we note that for our case, where 
c = c, equations (|4.6|), ( |4.7| ) and ( ^4.8| ) appear from equations (|4.3|), ( f4.4fl and Q4.5| ) by the 
interchangement a — > d. Hence the solution of the former three equations can immediately 
be read off from the solution of the latter three. 



The case ro = r\ = r 2 = r ^ 0: In this situation the equation ( [4.2[ ) gives that t 2 t\ = 1 
and the equations (|4.3|), (4.4) and (4^) give s\ = or s\ = t\ + t 2 . Thus we find the 



following four families of integrable models 
1. t 2 t\ = 1, si = s 2 = 0. 
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2. t 2 h = 1, si = 0, s 2 = + 

3. t 2 *i = 1, si=t± + t 2 , s 2 = 0. 

4. t 2 *i = 1, Si = s 2 = ti + t 2 . 

It is straightforward to check that these criteria for integrability do not depend on the 
choice of reference state. More precisely, the first three cases are related by Z% symmetry, 
i.e. they appear from one another when the labels 0, 1 and 2 are interchanged. The last 
is invariant under this symmetry. Thus, we have only two genuinely different classes of 
integrable models with r\ = r 2 = r% ^ 0. 

The case ro ^ 0, r\ = r 2 = 0: In this case it is obvious that the YBE's are always 
fulfilled since all the S-matrix elements are identical. 

The case ro = 0, r\ = r 2 = r / 0: Here the condition for integrability reads 

h = t 2 = si = s 2 . 
Finally, the model is obviously integrable for ro = n = r 2 = 0. 

5. R-matrices 

In this section we write down the .R-matrices corresponding to our integrable Hamiltonians 
with r\ = r 2 = r% / 0. These Hamiltonians all have their off-diagonal elements equal to 
one. They can therefore be characterized entirely in terms of their diagonal elements. We 
can choose representatives for the four different cases above for instance as follows 

1. H\ : (0, j, j,t,t + j, j,t,t,t + j), 

2. H 2 : (t + j, j, j,t, 0, j,t,t,t + j), 

3. H 3 : (t + j, j, j,t,t + j, j,t, t, 0), 

4. H 4 : (t + j,j,j,t,t + j,j,t,t,t + j), 

where the lists are lists of diagonal elements. Any member of a given class can be brought 
on the form above by addition of appropriate linear combinations of number operators. We 
can write the i?-matrix for the models above in a collective form as 

R oo(. u > v ) = A {u,v,s'), R\\(u,v) = A(u,v,s), R 2 2 l{u,v) = A(u,v, s), 

R\ 2 2 (u,v) = R 2 2 \(u,v) = R° 2 2 (u,v) = Rf (u,v) = R° \(u,v) = R\° (u,v) = C(u,v), 
R$(u,v) = R2 i(u,v) = R\\{u,v) = B(u,v,t), 

R%(u,v) = Rf (u,v) = R%(u,v) = B(u,v, ±), 
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where 

At \ su + uv + 1 

A(u,v,s) = — — , 5.1) 

sv + uv + 1 

1 + UV + jU + tv , 

B(u,v,t) = , (5.2) 

1 + UV + (j + t)v 

C(u,v)= , (5-3) 

1 + UV + (t + j)V 

and where the remaining .R-matrix elements vanish. The appropriate choices of s, s and s' 
for the four Hamiltonians above are 

1. #i : (s,s',s) = (i + ±,0,i + ±), 

2. H 2 : (s,s',s) = (t + i t+i 0), 

3. tf 3 : = (0,i+±,t+±), 

4. F 4 : (s,s',s) = (t + ±,t+±,f + ±). 

The R-matrix fulfills the necessary requirements, namely 

R(u,u)=P, (5.4) 
where P is the permutation operator and 

l + v 2 +(t+^)v^\ H = Pd u R(u,v)\ u=v . (5.5) 

Finally, it also satisfies the Yang Baxter relation 

R?tHu,v)RfiHu,w)R^Hv,w) = Rff(v,w)R^ 3 (u,w)R^ 2 (u,v). (5.6) 

As indicated by the notation above the i?-matrix has been constructed from pieces from 
the earlier determined S-matrix, cf. eqns. Q, (jig ), (pT|) , flO^ ), ( fug) and <$M) . 
In many situations one can perform a transformation from (u, v) to a new set of variables 
(A, u) such that R(\,v) = R(X — v). For instance, for the su(3) XXX spin chain, which 
belongs to case four and has t = 1, this transformation reads 

u = -——, v = — . (5.7) 

A + i v + i 

6. Comparison to known models 

Here we list a number of already known integrable models whose Hamiltonian can be 
written on the form (|2.1|). 



The first example of an integrable model of the form fl2.1| ) which comes to mind is 
the XXX su(3) spin chain. It is characterized by the elements of the Hamiltonian taking 
the values flg = H\\ = ff§ = 0, ffg = flg = H$ = flg = H™ = H% = 1, 
H 2i = H u = H oi = H w = H m = H 20 = - 1 or equivalents r = 1, s x = s 2 = 2, 
t\ = t 2 = 1. (Note that in the latter notation the angles are removed from the analysis.) 
It thus belongs to family number four. 
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In the same family we find the integrable deformations of this spin chain, describing 
the dilatation operator of the three scalar holomorphic sub-sector of the three parameter 
complex deformation of M = 4 SYM, given by (2.3) with r$ = n = r 2 = 1 |2C, [l8|]. They 
all have s\ = S2 = 2, t\ = t 2 = 1. 

Within our formalism we can of course investigate whether it is possible to achieve 
integr ability for a more general class of complex deformations of M = 4 SYM. The most 
interesting case is the case where all three r-variables are non- vanishing. First, we have 
seen that integr ability demands that ro = r\ = r 2 = r. In reference [^] it was argued that 
if ro = r\ = r 2 = r one furthermore needs that r = 1. This also follows from the results 
above. Namely, for the deformed model with ro = r\ = r 2 = r we have 

1 + r 2 , . 

si = s 2 = , (6.1) 

r 



2r 2 - 1 



(6.2) 



r 

2 — r 2 

h = , (6.3) 

r 

and we immediately see that according to the conditions for integrability presented above 
the model can only be integrable if r = 1 since only in this case t\t 2 = 1. 



In family four we also find the integrable su q (3) spin chain likewise studied in [20| for its 
possible connection to other deformations of J\f = 4 SYM. This model is characterized by 

ttOO — fill — FT22 _ n Wl2 — FfOl _ it02 _ i it12 _ it21 _ itIO _ trOl _ tT20 _ fj02 _ 
-"00 — -"11 — -"22 — u ' -"12 — -"01 — -"02 — "21 — "12 — "01 — "10 — "02 — "20 — 

r, H%1 = HIq = ilfo = 7,2 or equivalently s\ = S2 = h = r, t 2 = 1/r. 

In family three we find f. inst. the su(l|2) spin chain describing the su(l|2) sub-sector of 
undeformed N = 4 SYM and studied in |7|. Here one has H$ = ti,H\\ = 2, = 0, 

rr 12 _ rrl2 _ rr21 _ it 21 _ rrlO _ rrOl _ it20 _ it02 _ i irlO _ irOl _ it20 _ it02 _ i 
"12 — "21 — "12 — "21 — "10 — "01 — "20 — "02 _ L ' "01 — "10 — "02 _ "20 — L 

or equivalently r = 1, s\ = 0, s 2 = 2, t\ = t% = 1. The spin chain describing the analogue of 
the 5u(l|2) sub-sector in the three parameter complex deformation of J\f = 4 SYM theory 
is of course also included here. It has s\ = r + ^ — 2, S2 = r + ^, t\ = 2r — 7, £2 = 7 — r 



and is only integrable if r = 1 where it reduces to the model considered in [18]. 



We finally note some models known from studies of integrability in the context of con- 
densed matter systems, see for instance ^4j. Choosing the parameters as Hqq = ei cosh 7, 
H{\ = e 2 cosh 7 , EH = e 3 cosh 7 , 11% = = = -H% = -H$ = -H™ = sinh 7 
and H 2 l l = Hll = H$ = H§ = H$ = H$ = 1. Choosing e x = e 2 = e 3 = ±1 this is 



the anisotropic Perk-Schultz model [28] which then belongs to class 4. Setting 7 = in 



this model we find the su(3) Sutherland model [29]. Another example is the case when 
ei = — €2 = £3 = 1 which is the anisotropic supersymmetric t-J model [30], this model 
belongs to our class 2. 

7. Conclusion 

We have studied the most general three-state spin chain with nearest neighbour interaction 
and U(l) 3 symmetry. This model contains as a special case the spin chain describing the 
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three scalar holomorphic sub-sector of the three parameter complex deformation of Af = 4 
SYM, dual to type IIB string theory in the generalized Lunin-Maldacena backgrounds 
found by Frolov [15]. We have made use of the conceptually simple coordinate space Bethe 



ansatz invented by Bethe in 1931 22] and revived in connection with the study of M = 4 
SYM by Staudacher [23]. Solving the Yang Baxter relation for the iS-matrix we identified 
four classes of integrable models. Subsequently we wrote down an i?-matrix for the most 
interesting of these. Our findings show that each already known integrable model of the 
above type is nothing but one in a family of such models. We furthermore rule out the 
possibility that the complex three parameter deformation of the one-loop planar su(3) 
sector of M = 4 could be integrable under more general circumstances than those already 
identified. 
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A. The case of different dispersion relations 



The Yang Baxter relation (2.8) which describes the scattering between three excitatons 
has the well-known graphical interpretation shown in figure [|. A necessary condition for 
the relation to be fulfilled is that incoming and outgoing momenta are the same on the two 
sides of the equality sign. It is easy to find an example where this is not the case if r<i 7^ r\. 
For ro 7= one can f. inst. consider the situation shown in figure ||. Here the incoming 
particles are supposed to be of type 1, 1 and 2 with momenta qi, 52 and 53. The outgoing 
particles are assumed to be of type 2, 1 and 1 and their momenta can be found by the use 



of ( 3.25| ). It is easy to see that the outgoing momenta on the two sides of the equation do 



not match unless r 2 = T\. Since the integrability properties of the model can not depend 
on the choice of reference state the model can thus only be integrable if ro = r\ = r 2 . For 
ro = excitations of type 1 and 2 can not cross each but we can now in stead consider the 
diagram shown in figure ||. Again, we reach the conclusion that r 2 =T\. 

Let us finally comment on the case where the angles are present. In that case the 
analogue of ( 3.25|) reads: 

,1(72—71) p ipi+ip2 



e »p le -*(7i+72) T2 + rie " (A.l) 

i(-/2-'fl) p ipi+iP2 

(A.2) 



ip' 2 _ JP2j{l2+7i) ri + T2e e 



ri -I- r2 g2i(72-7i) e ipi+ip2 ' 

Setting r 2 = t\ we find 

Pi =Pi~ (72 + 7i)> P2 =P2 + (72 + 7l)- ( A -3) 
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Figure 1: Graphical representation of the Yang-Baxter relations 




1 



Figure 2: Factorized scattering requires this relation to hold. This is not possible unless the 
incoming and outgoing momenta are the same on both sides of the equality sign. Requiring that 
l'i = 41 implies that T2 =r\. 




i l 
Figure 3: In the case ro = this process determines that n = r%. 

This means that the outgoing momenta in for instance figure 2 are 
Qi = q" = <?l-(72 + 7l), 02 = <?2 = 92-(72 + 7l)> 93 = 93 = 93 + 2(72 + 7l)- (A.4) 

The outgoing momenta are hence the same on the two sides of the equation but the scat- 
tering is diffractive. However, if the system is integrable with the angles set to zero it is 
still integrable when the angles are introduced, cf. discussion on page ||. 
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